Elucidating the aryl hydrocarbon receptor antagonism from a chemical-structural perspective

The aryl hydrocarbon receptor (AhR) plays an important role in several biological processes such as reproduction, immunity and homoeostasis. However, little is known on the chemical-structural and physicochemical features that influence the activity of AhR antagonistic modulators. In the present report, in vitro AhR antagonistic activity evaluations, based on a chemical-activated luciferase gene expression (AhR-CALUX) bioassay, and an extensive literature review were performed with the aim of constructing a structurally diverse database of contaminants and potentially toxic chemicals. Subsequently, QSAR models based on Linear Discriminant Analysis and Logistic Regression, as well as two toxicophoric hypotheses were proposed to model the AhR antagonistic activity of the built dataset. The QSAR models were rigorously validated yielding satisfactory performance for all classification parameters. Likewise, the toxicophoric hypotheses were validated using a diverse set of 350 decoys, demonstrating adequate robustness and predictive power. Chemical interpretations of both the QSAR and toxicophoric models suggested that hydrophobic constraints, the presence of aromatic rings and electron-acceptor moieties are critical for the AhR antagonism. Therefore, it is hoped that the deductions obtained in the present study will contribute to elucidate further on the structural and physicochemical factors influencing the AhR antagonistic activity of chemical compounds

Extending Graph (Discrete) Derivative Descriptors to N-Tuple Atom-Relations

In the present manuscript, an extension of the previously defined Graph Derivative Indices (GDIs) is discussed. To achieve this objective, the concept of a hypermatrix, conceived from the calculation of the frequencies of triple and quadruple atom relations in a set of connected sub-graphs, is introduced. This set of subgraphs is generated following a predefined criterion, known as the event (S), being in this particular case the connectivity among atoms. The triple and quadruple relations frequency matrices serve as a basis for the computation of triple and quadruple discrete derivative indices, respectively. The GDIs are implemented in a computational program denominated DIVATI (acronym for DIscrete DeriVAtive Type Indices), a module of TOMOCOMD-CARDD program. Shannon‟s entropy-based variability analysis demonstrates that the GDIs show major variability than others indices used in QSAR/QSPR researches. In addition, it can be appreciated when the indices are extended over n-elements from the graph, its quality increases, principally when they are used in a combined way. QSPR modeling of the physicochemical properties Log P and Log K of the 2-furylethylenes derivatives reveals that the GDIs obtained using the tripleand quadruple matrix approaches yield superior performance to the duplex matrix approach. Moreover, the statistical parameters for models obtained with the GDI method are superior to those reported in the literature by using other methods. It can therefore be suggested that the GDI method, seem to be a promissory tool to reckon on in QSAR/QSPR studies, virtual screening of compound datasets and similarity/dissimilarity evaluations

Development of an in Silico Model of DPPH‚ Free Radical Scavenging Capacity: Prediction of Antioxidant Activity of Coumarin Type Compounds

A quantitative structure-activity relationship (QSAR) study of the 2,2-diphenyl-l-picrylhydrazyl (DPPH‚) radical scavenging ability of 1373 chemical compounds, using DRAGON molecular descriptors (MD) and the neural network technique, a technique based on the multilayer multilayer perceptron (MLP), was developed. The built model demonstrated a satisfactory performance for the training `R2 “ 0.713 ̆ and test set `Q2ext “ 0.654 ̆, respectively. To gain greater insight on the relevance of the MD contained in the MLP model, sensitivity and principal component analyses were performed. Moreover, structural and mechanistic interpretation was carried out to comprehend the relationship of the variables in the model with the modeled property. The constructed MLP model was employed to predict the radical scavenging ability for a group of coumarin-type compounds. Finally, in order to validate the model’s predictions, an in vitro assay for one of the compounds (4-hydroxycoumarin) was performed, showing a satisfactory proximity between the experimental and predicted pIC50 values

OncoOmics approaches to reveal essential genes in breast cancer: a panoramic view from pathogenesis to precision medicine

[Abstract] Breast cancer (BC) is the leading cause of cancer-related death among women and the most commonly diagnosed cancer worldwide. Although in recent years large-scale efforts have focused on identifying new therapeutic targets, a better understanding of BC molecular processes is required. Here we focused on elucidating the molecular hallmarks of BC heterogeneity and the oncogenic mutations involved in precision medicine that remains poorly defined. To fill this gap, we established an OncoOmics strategy that consists of analyzing genomic alterations, signaling pathways, protein-protein interactome network, protein expression, dependency maps in cell lines and patient-derived xenografts in 230 previously prioritized genes to reveal essential genes in breast cancer. As results, the OncoOmics BC essential genes were rationally filtered to 140. mRNA up-regulation was the most prevalent genomic alteration. The most altered signaling pathways were associated with basal-like and Her2-enriched molecular subtypes. RAC1, AKT1, CCND1, PIK3CA, ERBB2, CDH1, MAPK14, TP53, MAPK1, SRC, RAC3, BCL2, CTNNB1, EGFR, CDK2, GRB2, MED1 and GATA3 were essential genes in at least three OncoOmics approaches. Drugs with the highest amount of clinical trials in phases 3 and 4 were paclitaxel, docetaxel, trastuzumab, tamoxifen and doxorubicin. Lastly, we collected ~3,500 somatic and germline oncogenic variants associated with 50 essential genes, which in turn had therapeutic connectivity with 73 drugs. In conclusion, the OncoOmics strategy reveals essential genes capable of accelerating the development of targeted therapies for precision oncology.Instituto de Salud Carlos III; PI17/0182

Computational strategies for the discovery of biological functions of health foods, nutraceuticals and cosmeceuticals: a review

Scientific and consumer interest in healthy foods (also known as functional foods), nutraceuticals and cosmeceuticals has increased in the recent years, leading to an increased presence of these products in the market. However, the regulations across different countries that define the type of claims that may be made, and the degree of evidence required to support these claims, are rather inconsistent. Moreover, there is also controversy on the effectiveness and biological mode of action of many of these products, which should undergo an exhaustive approval process to guarantee the consumer rights. Computational approaches constitute invaluable tools to facilitate the discovery of bioactive molecules and provide biological plausibility on the mode of action of these products. Indeed, methodologies like QSAR, docking or molecular dynamics have been used in drug discovery protocols for decades and can now aid in the discovery of bioactive food components. Thanks to these approaches, it is possible to search for new functions in food constituents, which may be part of our daily diet, and help to prevent disorders like diabetes, hypercholesterolemia or obesity. In the present manuscript, computational studies applied to this field are reviewed to illustrate the potential of these approaches to guide the first screening steps and the mechanistic studies of nutraceutical, cosmeceutical and functional foods.We gratefully acknowledged the financial support of the Agencia Valenciana de la Investigación (AVI, https://innoavi.es/en/) by its program Innodocto (Reference number INNTAL32/19/002). SJB acknowledges the funding assistance from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 893810.Peer reviewe

Assessing the chemical-induced estrogenicity using in silico and in vitro methods

Multiple substances are considered endocrine disrupting chemicals (EDCs). However, there is a significant gap in the early prioritization of EDC’s effects. In this work, in silico and in vitro methods were used to model estrogenicity. Two Quantitative Structure-Activity Relationship (QSAR) models based on Logistic Regression and REPTree algorithms were built using a large and diverse database of estrogen receptor (ESR) agonism. A 10-fold external validation demonstrated their robustness and predictive capacity. Mechanistic interpretations of the molecular descriptors (C-026, nArOH,PW5, B06[Br-Br]) used for modelling suggested that the heteroatomic fragments, aromatic hydroxyls, and bromines, and the relative bond accessibility areas of molecules, are structural determinants in estrogenicity. As validation of the QSARs, ESR transactivity of thirteen persistent organic pollutants (POPs) and suspected EDCs was tested in vitro using the MMV-Luc cell line. A good correspondence between predictions and experimental bioassays demonstrated the value of the QSARs for prioritization of ESR agonist compounds

Triarylmethane class of compounds as potential modulators of the Aryl Hydrocarbon Receptor (AhR): synthesis and in vitro evaluation

editorial reviewedThe aryl hydrocarbon receptor (AhR) is a transcription factor and a crucial regulator of several functions in most cell types and tissues. The main outcome of AhR expression is the induction of cytochrome P450 family 1 (CYP1) of metabolizing enzymes, which has been recently suggested as a pharmacological target in endocrine homeostasis, cell-cycle regulation, and treatment of immune, inflammatory and oncological diseases, among others. Herein a novel set of triarylmethane (TAM) class of compounds is evaluated for possible modulatory activity of AhR expression. The biological evaluation of the AhR agonistic/antagonistic effects was performed by means of a cell-based bioassay for the set of TAMs. Similar to higher agonistic activity on AhR expression compared to the endogenous AhR agonist 6-formylindolo[3,2-b]carbazole (FICZ) was observed, and none of the TAMs tested showed any potential antagonism. Furthermore, an analysis of substituents influencing the activity of the bioisosters and predictions of the bioavailability parameters, the physicochemical properties and the drug-likeness of the novel TAMs were conducted.PROTECTion against Endocrine Disruptors; Detection, mixtures, health effects, risk assessment and communicatio

Predictive modeling of aryl hydrocarbon receptor (AhR) agonism

The aryl hydrocarbon receptor (AhR) plays a key role in the regulation of gene expression in metabolic machinery and detoxification systems. In the recent years, this receptor has attracted interest as a therapeutic target for immunological, oncogenic and inflammatory conditions. In the present report, in silico and in vitro approaches were combined to study the activation of the AhR. To this end, a large database of chemical compounds with known AhR agonistic activity was employed to build 5 classifiers based on the Adaboost (AdB), Gradient Boosting (GB), Random Forest (RF), Multilayer Perceptron (MLP) and Support Vector Machine (SVM) algorithms, respectively. The built classifiers were examined, following a 10-fold external validation procedure, demonstrating adequate robustness and predictivity. These models were integrated into a majority vote based ensemble, subsequently used to screen an in-house library of compounds from which 40 compounds were selected for prospective in vitro experimental validation. The general correspondence between the ensemble predictions and the in vitro results suggests that the constructed ensemble may be useful in predicting the AhR agonistic activity, both in a toxicological and pharmacological context. A preliminary structure-activity analysis of the evaluated compounds revealed that all structures bearing a benzothiazole moiety induced AhR expression while diverse activity profiles were exhibited by phenolic derivatives.PROTECTED - PROTECTion against Endocrine Disruptors; Detection, mixtures, health effects, risk assessment and communicatio

Atom Types Independent Molecular Mechanics Method for Predicting the Conformational Energy of Small Molecules

We previously implemented a well-known qualitative chemical principle into an accurate quantitative model computing relative potential energies of conformers. According to this principle, hyperconjugation strength correlates with electronegativity of donors and acceptors. While this earlier version of our model applies to σ bonds, lone pairs, disregarded in this earlier version, also have a major impact on the conformational preferences of molecules. Among the well-established principles used by organic chemists to rationalize some organic chemical behaviors are the anomeric effect, the alpha effect, basicity, and nucleophilicity. These effects are directly related to the presence of lone pairs. We report herein our effort to incorporate lone pairs into our model to extend its applicability domain to any saturated small molecules. The developed model H-TEQ 2 has been validated on a wide variety of molecules from polyaromatic molecules to carbohydrates and molecules with high heteroatoms/carbon ratios. Interestingly, this method, in contrast to common force field-based methods, does not rely on atom types and is virtually applicable to any organic molecules

DPPH• free radical scavenging activity of coumarin derivatives. In silico and in vitro approach

editorial reviewedThe interest of coumarins as antioxidant agents has attracted much attention in recent years. A quantitative structure-activity relationship (QSAR) study of the DPPH• (2,2-diphenyl-l-picrylhydrazyl) radical scavenging ability of chemical compounds, based on the 0-3D DRAGON molecular descriptors and an artificial neural networks (ANN) technique was developed. The built mathematical model showed a correlation coefficient for the training set (R2) = 0.71, an external correlation coefficient (2 ) = 0.65 and it was used to predict the antioxidant activity of 4-hydroxycoumarin derivatives. Besides, an experimental in vitro assay was developed for the reference compound of this group (4-hydroxycoumarin) and the results obtained confirmed the predictions made by the ANN